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ABSTRACT 

The kinematics of satellite galaxies reflect the masses of the extended dark matter 
haloes in which they orbit, and thus shed light on the mass-luminosity relation (MLR) 
of their corresponding central galaxies. In this paper we select a large sample of centrals 
and satellites from the Sloan Digital Sky Survey (SDSS) and measure the kinematics 
(velocity dispersions) of the satellite galaxies as a function of the r-band luminosity of 
the central galaxies. Using the analytical framework presented in Paper I, we use these 
data to infer both the mean and the scatter of the MLR of central galaxies, carefully 
taking account of selection effects and biases introduced by the stacking procedure. 
As expected, brighter centrals on average reside in more massive haloes. In addition, 
we find that the scatter in halo masses for centrals of a given luminosity, aiogM, 
also increases with increasing luminosity. As we demonstrate, this is consistent with 
fiogi, which reflects the scatter in the conditional probability function P{Lc\M), being 
independent of halo mass. Our analysis of the satellite kinematics yields a\ogL — 0.16± 
0.04, in excellent agreement with constraints from clustering and group catalogues, and 
with predictions from a semi-analytical model of galaxy formation. We thus conclude 
that the amount of stochasticity in galaxy formation, which is characterized by aiog l , 
is well constrained, is independent of halo mass, and is in good agreement with current 
models of galaxy formation. 

Key words: galaxies: haloes — galaxies: kinematics and dynamics — galaxies: fun- 
damental parameters — galaxies: structure — dark matter — methods: statistical 



1 INTRODUCTION 

According to the standard framework of galaxy formation, 
dark matter haloes form gravitational potential wells in 
which baryon s collapse, dissipate their energy and form star s 
and galaxies (| White fc Reesfl 19781 : iBlumenthal et al.lll984 ). 
The complex process of galaxy formation and evolution is 
believed to be governed by the mass of the dark matter halo 
in which it occurs. To understand the halo mass dependence 
of this process, it is important to statistically relate the ob- 
servable properties (e.g. luminosity) of galaxies to the masses 
of their dark matter haloes. 

The dark matter halo is bigger in extent than the vis- 
ible matter in a galaxy due to the inability of dark matter 
to lose energy via dissipative processes. Therefore, a tracer 
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population, distributed over large distances from the center 
of the halo is necessary to reliably measure the mass of a 
dark matter halo. Satellite galaxies satisfy this requirement 
and their kinematics reflect the depth of the potential well 
they orbit. However, a reliable measurement of the kine- 
matics requires a large number of satellites. Whereas a rel- 
atively large number of satellites can be detected in cluster- 
sized haloes, only a handful of satellites can be detected 
in less massive haloes. This leads to a large uncertainty on 
the kinematic measurement for less massive haloes. How- 
ever, under the assumption that central galaxies of similar 
properties reside in similar mass haloes, one can stack such 
central galaxies and use their satellites to obtain a statis- 
tical measure of the kinematics even in low mass haloes . 
lErickson, Gottesman & Huntei"l987'; 'Zar itskv et~allll993l : 
[Zaritsky & White 1994; Zaritskv ot al. 199'3)~^ 

During the past few years, data from large galaxy red- 
shift surveys such as the Sloan Digital Sky Survey (SDSS; 
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lYork et all bOOOD and the Two degree F ield Galaxy Red- 
shift Survey (2dFGRS: [Colless et al.ll200ll ') has become avail- 
able. The kinematics of satellite galaxies selected from 
these large datasets have been used to infer the mass- 
luminosity relationship (hereafter MLR) of central galaxies 
dMcKav et al.ll2002l: iBrainerd fc Specianll2003l; JPrada et af] 



I2OO3I : Ivan den Bosch et al.ll2004l : IConrov et all |2005|. l2007f ) 
and to study the density profiles of dar k matter haloes 
IjPrada et all I2OO3I : iKlvpin fc Pradal l2007l ). These studies 
consistently find that the velocity dispersion of satellites 
(and hence the mass of the halo) inc reases with the luminos- 
ity of the central galaxy. However, iNorberg. Frenk fc Cold 
l|2008l , hereafter N08) demonstrated a quantitative disagree- 
ment in the results of these studies and showed that this 
disagreement is largely due to subtle differences in the se- 
lection criteria used to identify central and satellite galaxies. 
Therefore, it is crucial to understand how selection effects 
bias the MLR of central galaxies inferred from satellite kine- 
matics. 

The MLR of central galaxies is specified in terms of 
the conditional probability P{M\Lc), which describes the 
probability for a central galaxy with luminosity Lc to reside 
in a halo of mass M. For a completely deterministic rela- 
tion between halo mass and central luminosity, P{M\Lc) = 
5[M — Mo), where 5 denotes the Dirac-delta function and 
Mo is a characteristic halo mass corresponding to centrals of 
luminosity Lc- The velocity dispersion, a{Lc), measured by 
stacking centrals with luminosity Lc, then translates into a 
mass Mo according to the scaling relation a'^ oc M. However, 
galaxy formation is a stochastic process and the distribution 
P{M\Lc) is expected to have non-zero scatter. If this scatter 
is appreciable then the stacking procedure results in com- 
bining the kinematics of haloes spanning a wide range in 
halo mass. This complicates the interpretation of the ve- 
locity dispersion. We addressed this issue in More, van den 
Bosch fc Cacciato (2008; hereafter Paper I), where we inves- 
tigated a method to measure both the mean and the scatter 
of the MLR of central galaxies using satellite kinematics. 
We outlined two different weighting schemes to measure the 
velocity dispersion of satellites, satellite-weighting and host- 
weighting, and showed that the mean and the scatter of the 
MLR can be inferred by modelling the velocity dispersion 
measurements in these two schemes simultaneously. In this 
paper, we apply this method to the spectroscopic galaxy 
catalogue from SDSS (Data Release 4) in order to deter- 
mine both the mean and the scatter of the MLR of central 
galaxies. 

The aim of this paper is twofold. First, we carry out a 
series of tests on a realistic mock galaxy catalogue to vali- 
date our method to infer the MLR of central galaxies from 
a redshift survey using satellite kinematics. In particular, 
we show that our central-satellite selection criterion and the 
method to measure the kinematics reliably recovers the true 
kinematics present in the mock catalogue. We also show that 
the mean and the scatter of the MLR inferred from the kine- 
matics match the corresponding true relations in the mock 
catalogue. Second, we repeat the analysis on galaxies from 
the SDSS. In particular, we show that both the mean and 
the scatter of the MLR increase with the luminosity of the 
central galaxy. This demonstrates that the scatter in halo 
masses is non-negligible and must be taken into account 
when interpreting measurements involving stacking. 



This paper is organized as follows. In Section [2] we de- 
scribe the construction of the mock catalogue that is used 
to test our method of analysis. In Section [3] we briefiy out- 
line the iterative selection criterion used to select centrals 
and satellites. In Section |4] we describe and test the method 
used to measure the kinematics of satellites as a function 
of the central luminosity. The inference of the MLR from 
the kinematics of satellite galaxies requires the knowledge of 
the number density distribution of satellites within a halo. 
In Section [5] we show that this distribution can also be in- 
ferred from the selected satellites. In Section |6] we describe 
our model to interpret the measured velocity dispersions and 
show that this model is able to recover the true mean and 
scatter of the MLR of central galaxies from the mock cata- 
logue. In Section [7] we apply our analysis method to a vol- 
ume limited sample from SDSS and discuss the results. We 
summarize our results in Section [8] 



2 MOCK CATALOGUES 

It is important to carefully identify central and satellite 
galaxies from a redshift survey in order to study the kine- 
matics of satellite galaxies. Furthermore, it is also important 
to reliably quantify the kinematics of the selected satellites 
as a function of central luminosity which in turn can yield 
the MLR of central galaxies. We monitor the performance 
of our method of analysis for each of these tasks using a re- 
alistic mock galaxy catalogue (MGC) which serves as a con- 
trol dataset. The halo occupation of galaxies in the MGC is 
known a priori, thereby allowing an accurate assessment of 
the level of contamination of the selected sample of centrals 
and satellites due to false identifications and also a com- 
parison between the kinematics recovered from the selected 
satellites and the actual kinematics present in the MGC. 

The two essential steps to construct a MGC are to ob- 
tain a distribution of dark matter haloes and to use a recipe 
to populate the dark matter haloes with galaxies. For the 
former purpose, we use a numerical simulation of dark mat- 
ter particles in a cosmological setup. For the latter, we use 
the conditional luminosity function (CLF) which describes 
the average number of galaxies with luminosities in the range 
L ± dI//2 that reside in a halo of mass M. 

A distribution of dark matter haloes is obtained from 
a A''-body simulation for a ACDM cosmology with param- 
eters, flm = 0.238, Ha = 0.762,0-8 = 0.75, n^ = 0.95, h = 
-ffo/100 kms^^ Mpc^^ = 0.73. The simulation consists of 
N = 512"^ particles within a cube of side Lbox = 300/i~'^Mpc 
with periodic boundary conditions. The particle mass is 1.33 
X 10^'^ h~^MQ. Dark matter haloes are identifi ed using the 
friends-of- friends algorithm (IDavis et al.l 119851 ) with a link- 
ing length of 0.2 times the mean inter-particle separation. 
Haloes obtaine d with this linking length have a m ean over- 
density of 180 (|Porciani. Dekel fc Hoffmanll2002l ). We con- 
sider only those haloes which have at least 20 particles or 
more. 

To populate the dark matter haloes with galaxies, we 
need to know the number and the luminosities of galaxies 
to be assigned to each halo. Furthermore, we also need to 
assign phase space coordi nates to each of these galaxies. We 
use the CLF described in lCacciato et al.l l|2008r ) for the first 
purpose. The CLF is a priori split into a contribution from 
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centrals and satellites, i.e. <I>(L|M) = $c(i|M) + <l>s(iiM). 
Here, ^c{L\M)AL denotes the conditional probability that 
a halo of mass M harbours a central galaxy of luminosity 
between L and L + dL, and $s(i|Af)dL denotes the average 
number of satellites of luminosity between L and L+dL. The 
parameters that des cribe the CLF are co nstrained using the 
luminosity function IjBlanton et al.ll2005r ) and th e luminosity 
depen dence of the correlation length of galaxies (|Wang et al.l 
l2007al ) in SDSS. 

Let us consider a halo of mass M. The luminosity of the 
central galaxy within this halo is sampled from the distribu- 
tion <E>c(L|Af ). The average number of satellites that have a 
luminosity greater than Lmin = 10^/i~^Lq and reside within 
haloes of mass M is given by 



{iV,at)(Af) 



$s(L|M)dL. 



(1) 



We assume Pois son statistics for th e occupation num- 
ber of satellites (iKrav tsov et all 12004 lYang et al.1 l2005al : 
lYang. Mo fc van den Bosch 200a) and assign A'sat galaxies 
to the halo where iVaat is drawn from 



P(iVsat|M)=exp(-M) 



A^sat! 



(2) 



with fj, = {Nsat){M). The luminosities of these satellite 
galaxies are drawn from the distribution $s(L|M). 

Phase space coordinates are assigned to the galaxies 
in the following manner. The central galaxy is assumed to 
reside at rest at the centre of the halo. Therefore, it has 
the same phase space coordinates as the parent dark matter 
halo. We assume that the halo is spherical and that the 
dark matter density distribution, p(r\M ), follows a universal 
profile (Navarro. Frenk fc Whitelll997l ') given by 



pWM)oc(i-)-^(l + i-)- 



(3) 



Here, rs denotes a scale radius which is specified in terms 
of the virial radius, rvir, using the concentration parame- 
ter, c(M) = ryjr/rs. We use the concentration-mass relation 
from lMaccio et al.l (2007) appropriately modified for our def- 
inition of the halo mass. The number density distribution of 
satellites, nsat(''|M), is given by 



nsat(r|M) oc (^^ 



n^ 



r 



7-3 



(4) 



where 7 represents the slope of the number density distri- 
bution of satellites, as r — » and 7?. is a free parameter. For 
populating the MGC, we adopt 7 = 7?. = 1 which implies 
that the satellites trace the dark matter density distribu- 
tion in an unbiased manner. The distribution, ns!it{r\M), is 
normalized such that 



{iVsat)(M)=47r 



nsat(r|M)r dr. 



(5) 



The radial coordinates of the satellite galaxies with respect 
to the center of the halo are sampled from the distribu- 
tion nsat(r|M). The satellite distribution around centrals is 
assumed to be spherically symmetric and random angular 
coordinates are assigned to the satellite galaxies. At the as- 
signed position for every satellite galaxy, one dimensional 
velocities are drawn from a Gaussian, 



f{Vj 



1 



/2¥a,at(r|M) 



exp 



■2'^Lir\M) 



(6) 



where Vj denotes the relative velocity of the satellite with 
respect to the central along axis j and (Tsat('"|Af) denotes 
the radial velocity dispersion at a distance r from the cen- 
tre of the halo. Here isotropy of orbits is assumed, i.e. the 
velocity dispersion along the j"^ axis, ct^ , equals crsat('"|A^)- 
The radial velocity dispersion within the halo is determined 
by solving the Jeans equation 



t{r\M)^ 



nsat(r|M) 



n.at(r'|M)|^(r'jAf)dr', (7) 



where ^(r'|Af) is the gravitational potential and the radial 
derivative represents the radial force given by 



-—{rM) = —^ 
or r^ 



p{r \M)r Ar . 



(8) 



We assume that the dark matter dominates the poten- 
tial. Using Eqs. JS} and (J4]) in Eqs. (O and JH]) gives 



2 

""sat 



(r\M) 



7^V(c) 



{— 



r 
7^^ 



3-7 



li{x)Ai 



(9) 



,^/^^ (a;/7^)^+2(l+x/7^)3-7' 
where 14ir is the circular velocity at rvir and 

^l{x)= j y{l + y)-^dy. (10) 

Jo 

The above expression for crsat(r'|Af) is used in the distribu- 
tion given by Eq. © to assign velocities to satellites. The 
entire procedure of assigning central and satellite galaxies is 
repeated for all the dark matter haloes within the simula- 
tion. 

Our aim is to construct a mock redshift survey 
that mimics the SDSS. Therefore, 2x2x2 identical galaxy- 
populated simulation boxes (which have periodic boundary 
conditions) are stacked together. A (RA, DEC) coordinate 
frame is defined with respect to a virtual observer at one of 
the corners of the stack. The apparent magnitude of each 
galaxy is computed according to its luminosity and distance 
from the observer. The line-of-sight (los) velocity of the 
galaxy is calculated by adding its peculiar velocity to the 
velocity of the cosmological fiow. A random velocity drawn 
from a Gaussian distribution with a dispersion of 35 kms~^ 
is further added along the los to account for the spectro- 
scopic redshift errors present in the SDSS. The redshift as 
seen by the virtual observer is then computed using the total 
velocity. We only consider galaxies with an observed red- 
shift z < 0.15 and an apparent magnitude brighter than 
17.77. This fiux limited catalogue is denoted henceforth by 
MOCKF and has 289,500 galaxies above an absolute lumi- 
nosity of lO^/i"^L0. MOCKF is used in Appendix [X] to 
investigate potential selection biases associated with the se- 
lection of central galaxies. In addition to MOCKF, we con- 
struct a volume limited sample, MOCKV, of galaxies that 
lie in the redshift range 0.02 ^ z ^ 0.072 and have luminosi- 
ties greater than 10^'^ H'^'Lq. It consists of 69,512 galaxies. 
In what follows, we use the volume limited sample MOCKV 
to validate our method for quantifying the kinematics (Sec- 
tion |4)|, to validate the method to infer the number density 
distribution of satellites (Section [5)l and finally to confirm 
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Table 1. Selection Criteria 



(AV)„ 

Figure 1. Schematic diagram of a selection criterion. Two coax- 
ial cylinders are defined around each galaxy (represented by a 
solid dot). The axis is along the los while the face of each cylin- 
der is parallel to the plane of the sky. 
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The parameters used to specify the inner and the outer cylin- 
der around a galaxy for the selection criterion used in this paper 
(ITER) and the one used in N08. The first row for ITER denotes 
the parameters used in the first iteration, while the second row 
denotes the parameters used in subsequent iterations. The veloc- 
ity dispersion, (7sat in units of 200kms~^ is denoted by (7200 and 
is used to scale the cylinders in every iteration. 



that the mean and scatter of the MLR can be reliably re- 
covered using the kinematics of satellites (Section [S|. 



3 SELECTION CRITERIA 

Large scale galaxy redshift surveys such as the SDSS allow 
the selection of a statistically significant sample of satellites. 
Since the observed galaxies cannot be a priori classified as 
centrals and satellites, it is important to use a selection cri- 
terion that can correctly identify central galaxies and the 
satellites which orbit around them. In this section, we de- 
scribe the selection criterion that we use to identify the cen- 
tral and satellite galaxies. 

A galaxy is identified to be a central if it is at least 
/h times brighter than every other galaxy within a cylindri- 
cal volume specified by _R < _Rh and |AV^| < (Al/)h (see 
Fig. [T|. Here, R is the physical separation from the candi- 
date central galaxy projected on the sky and AT/ is the los 
velocity difference. Around each of the identified centrals, 
satellites are those galaxies that are at least /s times fainter 
than their central galaxy and lie within a cylindrical volume 
specified hy R < Ra and |Ay| < (AV^)s. The identification 
of the central galaxies depends on the parameters Rh, {AV)h 
and /h, while the selection of satellites depends on the pa- 
rameters Rs, (A1/)s and fa. The values of these parameters 
also determine the level of contamination of the sample due 
to falsely identified centrals and falsely identified satellites 
(hereafter interlopers). The false identification of centrals 
can be minimized by choosing large values of i?h, (AT/)]! and 
/h so that the selected central is the dominant galaxy in a 
large volume. On the other hand, minimizing the interlopers 
requires small values of Rs and (AV)s. A large value of /s 
further guarantees that the selected satellites are small and 
do not dominate the kinematics of the halo (i.e. can safely be 
considered as test particles). Although stricter restrictions 
yield cleaner samples, they also reduce the sample size sig- 
nificantly. This makes the velocity dispersion measurements 
noisy. Thus, there is a tradeoff between the contamination 
level and the sample size. 

Most authors have chosen fixed values for the selec- 
tion criterion parameters, independent of the luminosity 
of the galaxy under consideration (McKay et al. 2002; 
Prada et al. 2003; Brainerd & Specian 2003; N08). Since 



brighter centrals on average reside in more extended haloes, 
Ivan den Bosch et al.l (|2004l ) advocated an aperture which 
scales with the virial radius of the halo around the galaxy. 
They used an iterative criterion which scales the cylindrical 
aperture based upon the estimate of the velocity dispersion 
around the central after every iteration. In this paper, we 
also use this iterative criterion to select centrals and satel- 
lites. The parameter set {Rh, (Al/)h, /h, Rs, (Ay)s, /a} that 
defines the inner and outer cylinders for the iterative crite- 
rion (ITER) is listed in Table [T] The first row lists the pa- 
rameters for the first iteration while the next row lists the 
scaling of these parameters in the subsequent iterations. In 
short, we proceed as follows: 

(i) Use fixed values of the aperture size to select centrals 
and satellites in the first iteration. 

(ii) Fit the velocity dispersion of the selected satellites as 
a function of the central galaxy luminosity, (Tsat(ic), with a 
simple functional form (see Section [4. l|l . 

(iii) Select new centrals and satellites by scaling the inner 
and the outer cylinders based on the estimate of the velocity 
dispersion. 

(iv) Repeat (ii) and (iii) until (Tsat(^c) has converged to 
the required accuracy. 

For step (iii), we a dopt the aperture scalings used in 
Ivan den Bosch et al.l (120041 ) . These aperture scalings were 
optimised to yield a large number of centrals and satellites, 
but at the same time reduce the interloper contamination. 
The values chosen for Rh and Rs approximately correspond 
to 2 and 0.375 times the virial radius, rvir- 



4 SATELLITE KINEMATICS 

In this section, we describe how to measure and model the 
velocity dispersion-luminosity relation, crsat(ic), using the 
satellites identified by the selection criterion. The relation 
o'sa.tiLc) can be measured either by binning the satellites by 
central galaxy luminosity or by using an unbinned estimator. 
We use the unbinned estimate after every iteration of the 
selection criterion to scale the selection aperture. However, 
to quantify the kinematics of the final sample of satellites, 
we use the binned estimator, for reasons which we describe 
further in the text. In the following subsections, we describe 
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Table 2. Parameters for the Selection Criterion 
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a 


b 


c 


MOCKV 


2.06 


0.45 


0.25 


MOCKF 


2.05 


0.50 


0.23 


SDSSV 


2.20 


0.38 


0.33 



The parameters used in Eq. JTTJ to define (T20o as a function 
of the luminosity of a galaxy in the final iteration for samples 
MOCKV, MOCKF and for the volume-limited sample from the 
SDSS, denoted by SDSSV. 



the unbinned and the binned estimators for (Tsat(ic) and 
finally an analytical model for the same. 



4.1 Unbinned Estimates 

We use a maximum likelihood method to estimate the rela- 
tion (Tsat(i'c) from the velocity information of the selected 
satellites after every iteration of the selection criterion. Let 
o"200 denote crsat(ic) in units of 200 kms~^ and Lio denote 
the lumin osity of the central galaxy i n units of 10^" h~'^\jQ. 
Following Ivan den Bosch et al.l l|2004 ). we parametrize (T200 
as, 



o-20o(logl/i 



a + 6 (log Lio) + c (log Li 



(11) 



Let /int denote the interloper fraction and assume that 
this fraction is independent of the luminosity of the cen- 
tral galaxy and AV^. The probability for a selected satellite 
to have a los velocity difference of A.V kms^^ with respect 
to the central is then given by 



P(AV) 



,/in 



2{^V).. 



where, o"eff 



+ 



+ ot 



/i. 



■ exp 



(AV)2 



2'^e'ff 



(12) 



11/2 



is the effective velocity disper- 



sion in the presence of the redshift errors and the factor 



27r(Toff erf 



{^V\ 



V2(Je 



(13) 



is such that the P(AV) is properly normalized to unity. 
In our attempt to mimic SDSS, we have added a Gaussian 
error of 35kms~^ to the velocity of each galaxy in the mock 
catalog. Therefore, the error on the velocity difference, AV, 
of the central and satellite galaxies is a^w ~ \/2 x 35kms~^ 
which adds in quadrature to (Jsat to yield (Jcff. 

We use Powell's direction set method to determine 
the parameters (a, b, c, /int) that maximize the likeli- 
hood 'C = ^. ln[P(Ay)]i, where the summation is over all 
the selected satellites. This yields a continuous estimate of 
o"sat(ic) without the need to bin the los velocity informa- 
tion of satellites according to the luminosity of the central 
galaxy. The parameter set (a, b, c) fitted in the last but one 
iteration determines the size of the apertures used to select 
the final sample of satellites. The values of these parame- 
ters for the samples investigated in this paper are listed in 
Table [3 
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Figure 2. Scatter plot of the velocity difference, AV, between 
the satellites and their centrals as a function of the central galaxy 
luminosity. The satellites were obtained by applying the iterative 
selection criterion to MOCKV. 



4.2 Binned Estimates 

We use a binned estimator to quantify the kinematics of the 
final sample of satellites. The binned estimator allows us 
to relax the simplistic assumption of /int being independent 
of Lc. More importantly, the binned estimator allows us, 
in a straightforward manner, to measure o-sat(ic) using two 
different weighting schemes - satellite-weighting and host- 
weighting. Most studies in the literature have used one of 
these two weighting schemes to infer the mean of the MLR. 
However, as demonstrated in Paper I, the mean of the MLR 
inferred from the velocity dispersion in any one of these 
two schemes is degenerate with the scatter in the MLR. 
This degeneracy can be broken by modelling the velocity 
dispersions in both schemes simultaneously. In what follows, 
we briefly explain these two weighting schemes in turn and 
then verify that the velocity dispersions in both schemes can 
be accurately recovered from the MGC. 

To measure the velocity dispersion of satellites in the 
satellite-weighting scheme, we obtain the distribution of ve- 
locities of the satellites, P(AV), with respect to their cen- 
trals for several bins of central galaxy luminosity (see Fig. 13]). 
Each bin has a width Alog[I/c] = 0.15. In this scheme, the 
centrals that have a larger number of satellites clearly con- 
tribute more to the P(A1/) distribution than those which 
have a smaller number of satellites. Therefore, the result- 
ing scatter in P(A1/) is a satellite-weighted average of the 
velocity dispersions around the stacked centrals (see Paper 
I for a detailed discussion). The dispersion obtained using 
this scheme is denoted henceforth by ctsw 

One has to undo the satellite- weighting described above 
in order to measure the host- weighted velocity dispersion. 
This can be accomplished by introducing a weight w — N^^ 
for each central-satellite pair while constructing the P{AV) 
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Figure 3. The satellite-weighted P{AV) distributions of satellites around centrals selected in several luminosity bins from MOCKV. 
The average log(Lc/ /i~^Lq) for each bin is indicate at the upper right corner of every panel. The (brown) dot-dashed line at the bottom 
of each distribution shows the contamination of the P{AV) distributions due to the interlopers. The (red) dotted lines indicate the single 
Gaussian fits, the (blue) dashed lines the double Gaussian fits. 



distribution jvan den Bosch et al.ll2004l : IConrov 6t"aLll2007l : 



iBecker et al.ll2007l ). Here, A'^ denotes the number of satellites 
selected around the central under consideration. Therefore, 
in this scheme each central receives a total weight of unity 
irrespective of the number of satellites it hosts. The scatter 
in this weighted P{AV) distribution is the host-weighted 
velocity dispersion and is denoted henceforth by (Jhw 

The procedure to obtain the scatter in the P{AV) dis- 
tributions is the same for both the satellite-weighted and 
the host-weighted case. This procedure must account for 
the interlopers and the redshift errors present in MOCKV. 
In what follows, we illustrate this procedure only for the 
satellite-weighted case. 

Fig.[2]shows the scatter plot of velocity difference AV of 
the selected satellites and the centrals as a function of the lu- 
minosity of the centrals. The satellite-weighted P{AV) dis- 
tributions of the satellites selected from MOCKV for several 
central luminosity bins are shown in Fig. [3] Dot-dashed lines 
show the contamination of the P{AV) distributions due to 
interlopers and are barely visible at the bottom of each dis- 



tribution. This confirms the claims in Ivan den Bosch et al.l 
(2004 ) that the iterative criterion yields a small fraction of 
interlopers with a weak dependence on Lc and that the in- 
terlopers can be modelled as a constant contribution to the 
velocity distribution independent of AV. 

A simple way to estimate the scatter of a P{AV) dis- 
tribution is to fit a Gaussian plus a constant model given 
by 



P(AV) ^ao + aiexp 






(14) 



Here, ao denotes the constant (with respect to AV) inter- 
loper background, ai is the normalization of the Gaussian 
and (Jeff is the effective dispersion in the presence of the 
redshift errors. 

The velocity dispersion obtained using a single Gaus- 
sian plus constant model fit can be systematically affected 
if the PjAV) d i stribu tion is intrinsically non-Gaussian. 
iDiaferio fc Gellerl (|l996l ) demonstrated that the velocity dis- 
tribution can be non- Gaussian partly due to mass mixing 
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Figure 4. Upper panels show the sateUite- weighted and the host- weighted velocity dispersions recovered from MOCKV. The (red) circles 
show values recovered from a single Gaussian fit while the (black) triangles show those from the double Gaussian fit. The solid line shows 
the variance of the true satellites and the dot-dashed line shows the analytical prediction using the halo occupation statistics of centrals 
from the MGC. The bottom panel plots the percentage deviation of the single and double Gaussian fits from the variance of the true 
satellites. 



(which is the result of stacking haloes of different mass) and 
partly due to the unrelaxed state of a halo. The second mo- 
ment of such a non- Gaussian distribution can be estimated 
with a double Gaussian plus a constant model (jBecker et al.l 
I2OO7I ) given by 



P{AV) =ao + aiexp 



2g\ 



4- a2exp 



\^y? 



2al 



The 


scatter 


(Tcff, in 


O-off 


aia'i 


+ a2(J2 


aicri 


+ a2(J2 



, in this case is such that 



-f (Je 



(15) 



(16) 



The dotted and the dashed lines in Fig. [3] show the single 
and double Gaussian fits to the P{AV) distributions respec- 
tively. 

Fig. 13] shows the velocity dispersions obtained from 
the satellite-weighted P(AV^) distributions in the upper left 
panel and those obtained from the host-weighted P{AV) 
distributions in the upper right panel. The (red) circles and 
the (black) triangles indicate the single and the double Gaus- 
sian fits respectively. Since the true satellites of centrals se- 
lected from MOCKV are known, they can be used to judge 
the goodness of the fits. The satellite- weighted and the host- 
weighted velocity dispersions of the true satellites (among 



the satellites selected using the iterative criterion) are ob- 
tained using 



2 






2 

ij 2 



(17) 



Here, N^ denotes the number of true centrals, Nj denotes 
the number of true satellites of the j* central and {AV)ij 
denotes the los velocity difference of the j^^ central with re- 
spect to its i"' satellite. The weight Wij = 1 for the satellite- 
weighted case and Wij — N~^ for the host-weighted case. 
The true velocity dispersions thus obtained are shown as 
solid curves in Fig. U) 

The bottom panels of Fig. [4] show the percentage devi- 
ation of both the single and the double Gaussian fits from 
the velocity dispersions of the true satellites. The single 
Gaussian fit (the dotted line) underestimates the dispersions 
systematically by about 5 — 10%. The double Gaussian fit 
(the solid line) on the other hand gives an unbiased esti- 
mate of both velocity dispersions. Therefore, in what fol- 
lows, we use the double Gaussian fit for measuring both the 
satellite-weighted and the host-weighted velocity dispersions 
(cf. Becker et al. 2007). 
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4.3 Analytical Estimates 

We now compare the velocity dispersions obtained from the 
satellite-weighted and the host-weighted schemes to their 
analytical expectation values. As detailed in Paper I, the 
satellite- weighted and the host-weighted velocity dispersions 
depend on the distribution of halo masses of central galaxies 
specified by P(M|Lc)- The analytical expressions describing 
the velocity dispersion in these two weighting schemes are 

k /an M KXIVI 

(18) 



2 .r ^ _ /o°° ^(A'/|ic)(Arsat)ap,A^(cTL)ap,MdM 



/q°° P(MjLc){iV,,t)ap,MdM 

2 /o°° P(M|Le) P({7V.at)ap,M) (aPap.Af AM 

"'"^ ^^ /o°°P(A/lL.)P({iV.at)ap,M)dM -^'^^ 

Here, the average number of satellites and the average ve- 
locity dispersion of satellites, within the aperture Rb in a 
halo of mass M, are denoted by {iVsat)ap,M and {crsat)ap,Af 
respectively. The number of satellites within the aperture, 
{A'"sat)ap,A/, is giveu by 



{iV.at) 



ap.JVI 



47r 



RAR 



nsat(r|M) 



rAr 



Vr 



7?2 



(20) 



Assuming the velocity dispersion of satellites to be isotropic, 
the velocity dispersion, (asat)ap,A/, can be expressed as the 
average of the radial velocity dispersion, (ysa.t{'''\^'^) (see 
Eq. [9} , over the aperture _Rs , 

{o-sat)ap,Af = TTT ^ / RAR 



{Asat)ap,A/ 

nsa.t{r\M) al, 



(r\M) 



rAr 



V~r 



R2 



(21) 



Note that, when measuring the host-weighted velocity dis- 
persions only satellites of those centrals that have at least 
one satellite within the search aperture are used. The frac- 
tion of such centrals is denoted by ■P((A'sat)ap,M) and is 
given by the probability that a halo of mass M, which on 
average hosts {A'sat)ap,A/ satellites within the aperture Ra, 
has A^sat ^ 1 within the aperture. Therefore, 



P({iVsat)ap,A./) 



P(Asat ^ 1) 

1 - P(Asat = 0) 

1 - exp[-{7V"sat)ap,M] 



(22) 



In the last equality we have assumed that the satellite occu- 
pation numbers (cf. Eq. [2]) follow Poisson s tatistics, which 
is su pported by numerical simulations (JKravtsov et al.l 
I2OO4I') and by resu l ts from group catalogs based on SDSS 
(JYang et ahlfeoOSal : IVang. Mo fc van den Boschll2008l ). The 
factor V{{Nsa,t)a.p,M) is not considered in the analytical es- 
timate in the satellite-weighting scheme as haloes with zero 
satellites, by definition, contribute zero weight. 

From the above analytical description, it is clear that 
the analytical estimates for the velocity dispersions require 
the knowledge of 

• the density distribution of dark matter haloes, p{r\M) 

• the number density distribution of satellites, naat(f |Af) 

• the halo occupation statistics of centrals, P{Al\Lc), and 
the halo occupation number of satellites, {Asat)(M). 

We assume that the density distribution of dark mat- 
ter haloes is given by Eq. ((3} and that the number density 



distribution of satellites is given by Eq. Q with 7 = 7?. = 1. 
The halo occupation statistics of centrals, P{M\Lc) is given 
by 



P(M|Lc) = 



$c(-Lc|M)n(M) 
J ^c{Lc\M)n{M)AM 



(23) 



where ^c{Lc\M) is the conditional luminosity function of 
central galaxies and n{M) is the halo mass function. The 
number of satellite galaxies in a halo of mass M is given 
by Eq. H}. We adopt the $c(I/c|M) and $s(L|Af) that 
were used in Section [2] to populate the MGC. With this 
input, we compute the analytical estimates for the veloc- 
ity dispersions of satellites as a function of luminosity using 
Eqs. (|18|l and (|19|) . The results thus obtained are shown 
as dot-dashed curves in the corresponding panels of Fig. U 
Overall the agreement with the velocity dispersions obtained 
from the satellites in the MGC is very good, except at inter- 
mediate luminosities where the analytical estimates are ~ 5 
percent higher than atruo. This indicates that the central 
galaxies selected from the MGC do not properly sample the 
full P{M\Lc). This can have two reasons: (i) a systematic 
problem with the criterion used to select central galaxies, 
or (ii) cosmic variance due to the finite volume probed by 
MOCKV. As we demonstrate in Appendix A our iterative 
criterion accurately samples the true P{AI\Lc), except for 
the fact that it misses the haloes of those centrals which have 
zero satellites. However, this sampling effect is accounted for 
in our analytical model via Eq. (|22p . In fact detailed tests 
show that the discrepancies between atiuc and our analytical 
estimates are entirely due to cosmic variance in the MGC. 

In Appendix [X] we also show that the strict selection 
criteria, that have been abundantly used in the literature, 
lead to a sample of central galaxies that is biased to reside 
in relatively low mass haloes. Consequently, the resulting 
MLR of central galaxies is similarly biased, and has to be 
interpreted with great care. 



5 NUMBER DENSITY DISTRIBUTION OF 
SATELLITES 

As described above, the number density distribution of satel- 
lite galaxies, nsat{r\M), is a necessary input to analytically 
compute the velocity dispersions. The projected number 
density distribution of satellites, E(7?|Lc), around centrals 
of a given luminosity, directly reflects the functional form of 
nsa.t{r\M). The distribution E(P|Lc) can be directly mea- 
sured by combining the satellites around centrals of a given 
luminosity, Lc, chosen by the selection criterion. However, 
it is necessary to first assess the impact of the interloper 
contamination on the measurement of E(7?|Lc), for which 
we again make use of the satellite sample selected from 
MOCKV. 

Fig. [5] shows, for the five brightest luminosity bins, 
the azimuthally averaged projected number density distribu- 
tions of the satellites selected from MOCKV. The errorbars 
refiect the Poisson noise on the number of satellites in each 
radial bin. The abrupt cutoff at large R is an artefact due to 
the parameter Rs in the selection criterion which describes 
the maximum projected radius within which satellites get 
selected. Note that, since Rs depends upon the luminosity 
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Figure 5. The projected number density distributions of the sateUites selected from MOCKV as a function of the projected radius in 
the brightest central luminosity bins. The errorbars assume Poisson statistics for the number of satellites in each bin. The (black) solid 
lines indicate the analytical predictions and assume that the satellite number density distribution follows the dark matter distribution in 
an unbiased manner, i.e. 7?. = 1 and 7 = 1 in eq. l|4]l. For comparison, the (red) dotted lines show the analytical predictions that assume 
7^ = 2 and 7 = 0. 



of central galaxies under consideration, this cutoff shifts to 
larger R with increasing central galaxy luminosity. 

The projected number density distribution of satellites 
around centrals stacked according to luminosity, E(_R|I/c), 
can be analytically expressed as. 



E(it:|Lc 



/ P(Af!Lc)E(7?|M)dM 

/ P(AflLc)P((Afsat)(M)ap,M) dM ' 



(24) 



Here, E(_R|M) is the projection of nsax{r\M) along the line- 
of-sight and is given by 



E(i?|M) 



■■' 7isat(r|M)2rdr- 



Vr 



B? 



(25) 



Using 7isat(''iAf) given by Eq. (|4]) with 7?. = 7 = 1 and the 
true P(M\Lc) present in the MGC, we analytically com- 
pute the expected number density distribution of satellites 
around centrals of a given luminosity. The solid lines in Fig. [5] 
show the results of this analytical expectation. The small dif- 
ferences between the measured and the analytically obtained 
distributions are due to the interlopers in the sample. How- 
ever, the differences become negligible in the brighter lumi- 
nosity bins. For comparison, the (red) dotted lines show the 
expected E(_R|Lc) for TZ = 2 and 7 = 0. This shows that the 
parameters TZ and 7, that characterize the number density 
distribution of satellites, can be inferred from the projected 
number density distributions of the selected satellites. 



6 MASS-LUMINOSITY RELATIONSHIP 

In the previous sections, using a mock catalog, we have 
demonstrated that the satellite-weighted velocity disper- 
sions, the host-weighted velocity dispersions and the pro- 
jected number density distributions of satellites around cen- 
trals of a given luminosity can be reliably measured starting 
from a volume limited redshift catalog of galaxies. Next, we 
attempt to infer the MLR of central galaxies from the ve- 
locity dispersions measured from MOCKV. The aim is to 
invert Eqs. (|18|l and (|19|l which describe the dependence of 
the velocity dispersions on the MLR of central galaxies. In 
addition to the velocity dispersions, we also measure the av- 
erage number of satellites per central of a given luminosity, 
{Nsa,t){Lc), and use this as a constraint. The dependence of 
{Nsat){Lc) on the MLR of central galaxies is given by 



(iVsat)(L,) = 



/^°°P(A/|Lc){iVsat)ap,AfdM 
/J" P(A./|Lc)P({iVsat)ap,M)dM' 



(26) 



In this section, we first describe the model we use to infer 
the mean and the scatter of the MLR from the observables 
o"sw, o"hw and (A^sat). Next, we use this model to infer the 
mean and the scatter of the MLR in the MGC and compare 
it to the true relations present in the MGC. 



6.1 The Model 

As mentioned earlier, the analytical computation of ctsw, cthw 
and (A'aat) requires the knowledge of the density distribu- 
tion of dark matter haloes, the number density distribution 
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Figure 6. Crosses with errorbars in the upper panels denote the data used to constrain the Monte Carlo Markov Chain; the satellite- 
weighted velocity dispersions in panel (a), the host- weighted velocity dispersions in panel (b) and the average number of satellites 
per central in panel (c). The relations recovered from the MCMC analysis are shown in the bottom panels; {log Lc){M) in panel (d), 
(logM)(Lc) in panel (e) and criog jv/(Lc) in panel (f). In each panel, the blue and purple colours denote the 68% and the 95% confidence 
levels. The solid lines in the lower panels denote the true relations present in MOCKV. 



Table 3. MOCKV: Parameters recovered from the MCMC 



Parameter Input 16% 50% 84% 



log(-Lo) 
log(Mi) 

72 



9.93 
11.04 
0.25 
0.14 



9.64 
10.48 
0.18 
0.13 



10.01 10.32 

11.28 11.69 

0.26 0.32 

0.15 0.17 



The input parameters that describe P{Lc\M) are compared to 
the 16"^, 50**^ and the 84**^ percentiles of the corresponding dis- 
tributions of parameters obtained from the MCMC. 



of satellites and the halo occupation statistics of centrals and 
satellites. We assume that the density distribution of dark 
matter haloes follows the NFW profile given by Eq. Q . For 
the number density distribution of satellites within a halo 
of mass M, nsa.t{r\M), we use Eq. (|4)| with 7?. = 7 = 1. As 
shown in Section [S] the projected number density distribu- 
tions of satellites selected from MOCKV is consistent with 
this analytical expression. Next, we describe our model for 
the halo occupation statistics of the centrals, specified by 
P(M|Lc), and the satellites, specified by {Nss,t}{M). 

The distribution P{M\Lc) is related to the complemen- 
tary distribution, P{Lc\M), by Bayes' theorem 



P(M|Lc) = 



n{M)P{Lc\M) 
Jn{M)P{Lc\M)dM' 



(27) 



where n{M) i s the halo mass function. We follow 
ICacciato et al.l (|2008l ) and parametrize the distribution 
P{Lc\M) as a lognormal in Lc, 



P(Lc|M)dLc 



log(e) 



27rcriog L 



■ exp 



(log[Lc/Lc])' 



^<L 



dLc 



(28) 



Here logLc(M) denotes the mean of the lognormal distribu- 
tion and (JiogL is the scatter in this distribution. We use four 
parameters to specify the relation Lc{M): a low mass end 
slope, 71, a high mass end slope, 72, a characteristic mass 
scale, M\, and a normalisation, Lq; 



Lc — Ln 



{M/Mi 



[l-h(M/Mi)]Ti-T2' 



(29) 



We assume the scatter uiogL to be independent of mass. 
We do not explore the faint end slope, 71, in our anal- 
ysis as the velocity dispersions at the faint end are very 
uncertain due to low number statistics. Instead, we keep 
it fixed at 3.273, which is the value obtained from the 
analysis of the abu ndance and clustering of galaxies (see 
ICacciato et al.ll2008lV This parametriz a tion i s motivated by 
results of lYang. Mo fc van den BoschI (|2008r ) who measure 
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the conditional lumin osity function from the SDSS group 
catalogue described in lYang et alj (|2007l ). 

We model the satellite occupation number, {Ns!it}{M), 
as a power law distribution, given by 



{iV,at)(A/)=iV. 



M 



1012 /i-1Mq 



(30) 



which adds two more parameters {Ns,a). Thus, in total, 
our model has six free parameters (aiog l,Lo, Mi ,y2,Ns,a). 
Given these parameters and the radial number density dis- 
tribution of satellites (specified by 7?. and 7), the velocity 
dispersions (Tsw(Lc) and (Thw(ic) as well as the number of 
satellites per central, {Nsat){Lc) in an aperture of a given 
size can be computed using Eqs. (|18|l . (|19p and (|26p and 
compared to the measured values. Crosses with errorbars in 
Panels (a), (b) and (c) of Fig. [6] show CTsw, cri™ and (A^sat) 
as a function of the luminosity of the central obtained from 
MOCKV. We use these measurements to constrain the six 
free parameters of our model. 



6.2 Monte-Carlo Markov Chain 

To determine the posterior probability distributions of the 
6 free parameters in our model, we use the Monte-Carlo 
Markov Chain (hereafter MCMC) technique. The MCMC 
is a chain of models, each with 6 parameters. At any point 
in the chain, a trial model is generated with the 6 free pa- 
rameters drawn from 6 independent Gaussian distributions 
which are centred on the current values of the correspond- 
ing parameters. The chi-squared statistic, Xtry, for this trial 
model, is calculated using 



2 

Xtry = 

with 



2 

Asw 



Ahw 



2 

Xns 



2.2.2 
Xsw ~r Ahw ~r Ans 



E 

1=1 

10 

E 

1=1 

10 
E 



Crsw(ic[J]) " O-sw(icW) 



Aas„{Lc[i\) 




^hv,{Lc[i]) -o-hw(ic[i]) 


2 


Aahw(ic[i]) 




{iVBat)(icW)-iV,,t(icW) 



AiV,at(I'c[i]) 



(31) 

(32) 

(33) 

(34) 
(35) 



Here X denotes the observational constraint X and AX 
its corresponding error. The trial step is accepted with a 
probability given by 



accept 



1.0, 



exp[ 



-ix' 



try 



Xcur 



)/2], 



ry 



$: 



Xcur 
-"-Xtry ^ Xcur 



ifx? 



(36) 



where Xcur denotes the x^ for the current model in the chain. 
We initialize the chain from a random position in the 
parameter space and discard the first lO** models allowing 
the chain to sample from a more probable part of the distri- 
bution. This is called the burn-in period for the chain. We 
proceed and construct a chain of models consisting of 10 mil- 
lion models. We thin this chain by a factor of 10* to remove 
the correlations between neighbouring models. This leaves 
us with a chain of 1000 independent models that sample the 



posterior distribution. We use this chain of models to esti- 
mate the confidence levels on the parameters and relations 
of interest. 

In Table El we compare the 16*, 50"" and the 84*'' 
percentiles of the distributions of parameters, which char- 
acterize P(Lc|M), obtained from the MCMC with the 
corresponding true values of these parameters present in 
MOCKV. The true parameter values have been recovered 
within the 68% confidence intervals. The 68 and 95% con- 
fidence levels in panels (a) , (b) and (c) of Fig. |6] show that 
the models from the MCMC accurately fit the velocity dis- 
persions, (Jsw and (Thw, as well as the average number of 
satellites per central, (Nsat) as a function of central galaxy 
luminosity. The confidence levels for the average luminos- 
ity of the centrals as a function of the halo mass, Lc{M), 
are shown in panel (d). The confidence levels on the mean, 
(log A4){Lc), and the scatter, (JiogM{Lc), of the MLR of cen- 
tral galaxies are shown in panels (e) and (f ) respectively. The 
solid lines in the lower panels show the corresponding true 
relations present in MOCKV. Clearly, our method is able to 
accurately recover the true MLR. 

This completes our tests with the MGC. Employing a 
variety of tests on a realistic MGC, we have established 
a proof-of-concept that, starting from a redshift survey of 
galaxies, one can reliably select central and satellite galax- 
ies, quantify the kinematics of the selected satellites around 
central galaxies and use this information to infer an unbi- 
ased estimate of the mean and the scatter of the MLR of 
central galaxies. 



7 RESULTS FROM SLOAN DIGITAL SKY 
SURVEY ANALYSIS 

We now apply the method tested in the previous sections to 
data from the Sloan Digital Sky Survey (SDSS; lYork et al.l 
|2Q0(J ) which is a joint five-passband {u,g,r,i and z) imag- 
ing and medium resolution {R ~ 1800) spectroscopic survey. 
More specifically, we use the New York University Value 



Added Galaxy Catalogue (jBlanton et al.l 12005), which is 



based upon SDSS Data Release 4 ( Adelman-McCarthv et al.l 
[200a) but includes a set of significant improvements over 
the original pipelines. From this catalog, we select all galax- 
ies in the main galaxy sample with redshifts in the range 
0.02 ^ z ^ 0.072 and with a redshift completeness limit C > 
0.8. We construct a volume limited sample of galaxies that 
have a r-band luminosity (k-corrected to redshift 0.1) above 
^min = 10®'^ h~^LQ. This sample is henceforth denoted by 
SDSSV. It consists of 57, 593 galaxies. The SDSSV data is 
analysed assuming the cosmol ogical parameters fro m the 3 
year data release of WMAP IjSpergel et all l2007h . Qm = 
0.238, Oa = 0.762, h = Ho/lOO kms"^ Mpc'^ = 0.734, 
the spectral index of initial density fiuctuations ris — 0.951 
and normalization as = 0.744. 

The iterative criterion (ITER) outlined in Section [3] is 
applied to select centrals and their satellites from SDSSV. 
The parameters (a, b, c) in Eq. (Ilip that define the aperture 
used in the final iteration of the central-satellite selection are 
listed in the last row of Table (2] The total number of central 
galaxies that host at least one satellite is 3863. The number 
of satellite galaxies selected is 6101. 

Fig. [7] shows the scatter plot of the velocity difference. 
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Figure 7. Scatter plot of the velocity diflference, AV, between 
the satellites and their central galaxies as a function of the central 
galaxy luminosity. The central galaxies and satellites are selected 
from SDSSV using the iterative selection criterion described in 
Section |3] 
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Figure 8. The satellite-weighted (red triangles) and the host- 
weighted (blue squares) velocity dispersions obtained from cen- 
trals and satellites selected from SDSSV. 



Table 4. SDSSV: Velocity Dispersion Measurements 



log(Lc) 


0"sw 


Actsw 


o"hw 


Ao-hw 


h-^LQ 


kms~^ 


kms"-"^ 


kms"-"^ 


kms-i 



9.61 


108 


20 


107 


20 


9.73 


148 


26 


146 


25 


9.88 


159 


19 


155 


18 


10.03 


162 


14 


159 


12 


10.17 


214 


43 


203 


20 


10.31 


254 


11 


220 


11 


10.45 


272 


15 


247 


13 


10.59 


412 


35 


287 


39 


10.72 


470 


28 


378 


73 


10.87 


650 


54 


574 


254 



The velocity dispersion measurements in the satellite-weighted 
and host-weighted schemes together with the associated errors. 



AV, between the selected satellites and the corresponding 
centrals as a function of the central luminosity. The scatter 
in the velocities of satellites with respect to their centrals 
clearly increases with central galaxy luminosity. To quantify 
this scatter, we obtain the P{AV) distributions in both the 
satellite-weighting and the host-weighting schemes by com- 
bining the velocity differences, A V, of satellites within lumi- 
nosity bins of uniform width A log[I/c] = 0.15. The satellite- 
weighted and the host-weighted velocity dispersions are es- 
timated from these distributions by fitting a double Gaus- 
sian plus constant model as described in Section 14.21 Fig. [8] 
shows these dispersions as a function of central luminos- 
ity. The values of ctsw, criiw and their associated errors are 
listed in Table U Both the satellite- weighted and the host- 
weighted velocity dispersions increase with the luminosity of 
the central galaxy. Note that the satellite-weighted velocity 
dispersions are systematically higher than the host- weighted 
velocity dispersions. As is evident from Eqs. psp and (|19p . 
this is a sufficient condition to indicate the presence of scat- 
ter in the MLR of central galaxies (see Paper I for a detailed 
discussion) . 

The model to infer the MLR of central galaxies from the 
kinematics requires the radial number density distribution 
of satellites, risat (r) , as an input. For inferring the MLR from 
the kinematics measured from MOCKV, we used a model 
of risat (r) that follows the density distribution of the dark 
matter in an unbiased manner, i.e., 7 = 7?. = 1 in Eq. Q. 
However with SDSS, it is not clear what functional form of 
nsa.t{r) should be used. In fact, various studies have shown 
that the satellite galaxie s are spatially ant i biased with re- 
spect to the dark matter (|Yang et al.ll2005bl : IChenll2007al lbl). 
Rather than including 7 and 7?. as free parameters in our 
model, we seek to constrain these parameters using the ob- 
servable E(i?|Lc). Fig.l^shows the projected number density 
distributions of the selected satellites for the five brightest 
luminosity bins. As can be seen from Eq. (|24p . predicting 
E(_R|Lc) requires the knowledge of P{M\Lc), which is the 
principle goal of our study. It further also requires the knowl- 
edge of {Ns!it){M). Both these quantities are unknown. We 
proceed as follows . We use the P(M\Lc] and (Afsat)(M) from 
the CLF model of lCacciato et al.l (|2008l ) which was also used 
to populate the mock catalogue. We explore two different 



models for naa.t{r), one with 7?. = 7 = 1, where the number 
density distribution of satellites follows the dark matter den- 
sity distribution, and the other with TZ = 2 and 7 = 0, where 
the number density distribution of satellites is spatially an- 
tibiased with respect to the dark matter distribution. The 
former model is shown as a (black) solid line while the latter 
with a (red) dotted line in Fig. [9] Clearly, the latter model 
is favoured by the data. Therefore, we use TZ = 2 and 7 = 
to specify nsat (r) for the analysis of the velocity dispersions 
to infer the MLR of central galaxies. 

Next we use the parametric model described in Sec- 
tion l6.1l and constrain it using the measured velocity disper- 
sions, (Taw and (Thw, and the number of satellites per central, 
(A'sat) as a function of the luminosity of centrals. This allows 
us to determine both the mean and the scatter of the MLR. 
Using a MCMC, we recover the mean relations (log M){Lc), 
mosM{Lc) and (logLc)(M). The 16*, 50*'^ and the 84*'^ per- 
centiles of the distributions of the parameter values that de- 
scribe the P{M\Lc) distribution for the central galaxies are 
listed in Table [S] Fig. [10] is the equivalent of Fig. [6] but for 
SDSSV. The upper panels show the data used to constrain 
the parameters and the bottom panels show the confidence 
levels on the inferred mean and scatter of the MLR. The 
values of (logM)(Lc) and (TiogAf(ic) together with their 1- 
(j errors are listed in Table [B] 

Clearly the average masses of dark matter haloes in- 
crease with central galaxy luminosity, as expected. Interest- 
ingly, the scatter in halo masses also increases systematically 
with the luminosity of the central galaxy. At the bright end, 
this scatter is roughly half a dex. Therefore, stacking central 
galaxies by luminosity amounts to stacking haloes that cover 
a wide range in masses. This justifies the need to account 
for this scatter in the analysis of the satellite kinematics. 
Neglecting this scatter leads to an overestimate of the halo 
mass at a given central luminosity. Most previous studies 
dealing with satellite kinematics have neglected this scatter 
which has resulted in a biased estimate of the halo mass- 
luminosity relationship. As we show in Appendix \^ their 
use of strict selection criteria to identify the centrals and 
satellites have further biased their estimate of the MLR of 
central galaxies. 



7.1 Comparisons with Independent Studies 

In a recent study, ICacciato et al.l (|2008l ') have constrained 
the CLF using the abundance and clustering of galaxies in 
SDSS. They have shown that this CLF is also able to re- 
produce the galaxy-galaxy lensing signal and is further con- 
sistent with the MLR obtained fro m a SDSS group catalog 
I Yang. Mo fc van den Boschll2008l ). It is interesting to com- 
pare the results of their study with the results obtained here 
from satellite kinematics. The solid lines in the bottom pan- 
els of Fig. [To] show the mean and the scatt er of the MLR 
obtain ed using the best-fit CLF model from ICacciato et al.l 
112003) . Clearly, these results are in excellent agreement with 
the results obtained here using the kinematics of satellite 
galaxies. Amongst others this provides further support that 
the halo mass as signment i n the SPSS group catalog of 
JYang et all (|2007l ) is rehable (|Wang et al.ll2007ibl ). 

Since the mean and scatter of the MLR reflect the 
physics, and in particular the stochasticity, of galaxy forma- 
tion, it is interesting to compare the results obtained here 
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Figure 11. The posterior distribution of the parameter cri^g^ 
as obtained from the MCMC analysis of the satellite velocity 
dispersions. The l-c constraints on the parameter (JiogL ob- 
tained from other independent methods are shown as shaded re- 
gions. Region GC ind icates the SDSS group catalog result by 
lYang. Mo fc van den B osch (2008)), region CLF indicates the re- 
sult obtained by ICooravl (120061 ) with an independent CLF anal- 
ysis and region SAM sh ows our measuremen t of criog j^ from the 
semi-analytical model of lCroton et al.l 1120061 1. 



to predictions from semi-analytical models (S AM) of galaxy 
forma tion. To that extent we use the SAM of ICroton et al.l 
(J2006I ). which has been shown to match the observed prop- 
erties of the local galaxy population with reasonable accu- 
racjo. Using a volume limited sample of galaxies selected 
from the SAM with the same luminosity and redshift cuts 
as SDSSV, we measure the mean and the scatter of the 
distributions of halo masses for central galaxies in several 
bins of r-band luminosity. The results are shown in panels 
(e) and (f ) of Fig. [10] as open squares. The agreement with 
our constraints from the satellite kinematics is remarkably 
good. It is both interesting and encouraging that a semi- 
analytical model, which uses simple, physically motivated 
recipes to model the complicated baryonic physics associ- 
ated with galaxy formation, is able to reproduce not only 
the mean of the MLR of central galaxies but also the cor- 
rect amount of stochasticity in this relation. 

In our model, the stochasticity of galaxy formation 
is best described by the parameter aiogi, which indicates 
the amount of scatter in the luminosity of central galaxies 
given the mass of a halo, i.e. the scatter in the distribu- 
tion P{Lc\M). The histogram in Fig. llll shows the posterior 
probability of aiog l, obtained from our Monte Carlo Markov 
Chain, which yields that aiogL = 0.16±0.04 (68% confidence 
levels) . Note that we have made the assumption that aiog l is 
independent of halo mass. The same assumption was made 



^ Note that ICroton et al.l 1120061 ) adopted a slightly different cos- 
mology than the one used in our data analysis which can have a 
small impact on the MLR. 
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Figure 9. The projected number density distributions of satellites around centrals selected from SDSSV for the five bright luminosity bins 
(average log(Lc/ /i^^Lq) in right corner). The (black) solid curves indicate the expected distributions if the number density distribution 
of satellites follows the dark matter density, i.e. 7^ = 7 = 1 in Eq. (|4]l. The (red) dotted curves in turn indicate the expected distributions 
for a model in which the satellite galaxies are a factor 2 less concentrated than dark matter and have a central core in the number density 
distribution i.e. 7?. = 2, and 7 = 0. 



by ICooravl ^200(1), who obtained that ai^gL = 0.17+o;°i 
using the luminosity f unction and clusterin g properties of 
SDSS galaxies (see also ICacciato et al.l 20081) . Using a large 
SPSS galaxy group catalogue, lYang. Mo fc van den BoschI 
(I2OO3) obtained direct estimates of the scatter in P(Lc|M), 
and found that criog_t = 0.13 ± 0.03 with no obvious depen- 
dence on h alo mass. Fina Uy, we also determined criogL in 
the SAM of ICroton et al.l ([2006): using several bins in halo 
mass covering the range 10^° /I'^M© < M ^ 10^^ h"^MQ, 
we find that aiogi = 0.17 ± 0.02, once again with virtually 
no dependence on halo mass. All these results are summa- 
rized in Fig. 111! Not only do they support our assumption 
that aiogL is independent of halo mass, they also are in re- 
markable quantitative agreement with each other. 

Finally we note that a constant scatter, (JiogL, in the 
distribution P{Lc\M) leads to a scatter, aiog ai(Lc), in the 
distribution P{M\Lc) that increases systematically with the 
luminosity of the central galaxy. This arises from the fact 
that the relation {logLc)(M) (Panel (d) in Fig. (TO]) is shal- 
lower at the high mass end compared to the low mass end 
(see Paper I). 



8 SUMMARY 

The kinematics of satellite galaxies have been widely used 
to statistically relate the mean halo masses of central galax- 
ies to their luminosities (Zaritsky et al. 1993; Zaritsky & 
White 1994; Zaritsky et al. 1997; McKay et al. 2002; Prada 
et al. 2003; Brainerd & Specian 2003; Norberg et al. 2008). 



Table 5. SDSSV: Parameters recovered from the MCMC 



Parameter 16% 50% 84% 



log(i()) 


9.69 


10.05 


10.33 


log(Mi) 


10.85 


11.74 


12.01 


72 


0.19 


0.28 


0.35 


0"logL 


0.12 


0.16 


0.19 



The 16*1^, 50"^ and the 84"^ percentiles of the distributions of 
parameters that describe the relation Lc{M) obtained from the 
MCMC analysis of the velocity dispersions obtained from SDSSV. 

These studies use strict criteria to identify central and satel- 
lite galaxies that reside preferentially in iso lated environ- 
ments. Following v an den Bo s ch et al.l (|200J), we applied a 
relaxed but adaptive selection criterion to a volume limited 
sample from SDSS to identify centrals and their satellites, 
not only in isolated environments but also in massive groups 
and clusters which allowed us to study the kinematics of 
satellites over a wide range of central galaxy luminosity. 
We inferred both the mean and the scatter of the mass- 
luminosity relationship of central galaxies from the kinemat- 
ics of satellite galaxies. 

Using a realistic mock catalogue, we thoroughly tested 
the analysis method at every step. We tested the perfor- 
mance of our criterion to identify central and satellite galax- 
ies and our method to measure the kinematics of the selected 
satellites. We showed that the kinematics recovered from the 
selected satellites are a fair representation of the true kine- 
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Figure 10. Crosses with errorbars in the upper panels show the data points used to constrain the MCMC; the satellite-weighted velocity 
dispersions in panel (a), the host- weighted velocity dispersions in panel (b) and the mean number of satellites per central as a function 
of luminosity in panel (c), all measured by using the satellite sample selected from SDSSV. The blue and purple bands represent the 
68% and 95% confidence regions respectively. The bottom panels show the relations inferred from the MCMC; the average log(Lc) is 
in panel (d) , and the mean and the scatter in the MLR of central galaxies in panels (e) and (f ) respectively. The relations obtained by 
using the halo occupation distribution fro m MOCKV are show n using solid lines. The squares in panels (e) and (f) indicate the values 
obtained from the semianalytical model of lCroton et al.l 112009) . 



Table 6. SDSSV: MLR of Central Galaxies 



log(Lc) 


(logM) 


A(logM) 


0"logM 


Ao"log M 


fc-2L0 


h-^MQ 


h-^MQ 


h-^M0 


h-^Me 


9.61 


12.06 


0.35 


0.12 


0.06 


9.73 


12.16 


0.32 


0.13 


0.07 


9.88 


12.28 


0.29 


0.15 


0.08 


10.03 


12.44 


0.26 


0.18 


0.10 


10.17 


12.60 


0.23 


0.22 


0.10 


10.31 


12.80 


0.21 


0.26 


0.11 


10.45 


13.01 


0.19 


0.30 


0.10 


10.59 


13.24 


0.19 


0.34 


0.09 


10.72 


13.47 


0.21 


0.36 


0.08 


10.87 


13.74 


0.23 


0.38 


0.07 



The mean and scatter of the halo masses as a function of the 
central galaxy luminosity inferred from the MCMC analysis. The 
errors on each of the inferred quantities correspond to the 68% 
confidence levels. 

matics of satellite galaxies present in the mock catalogue. 
We presented an analytical model that properly accounts 
for the selection biases and showed that the predictions of 



this analytical model are in good agreement with the mea- 
sured kinematics of the selected satellites. 

In Paper I, we have shown that the velocity disper- 
sion of satellites can be measured using two different weight- 
ing schemes: satellite- weighting and host- weighting. We have 
demonstrated a degeneracy between the mean and the scat- 
ter of the MLR obtained from either the satellite-weighted 
or the host-weighted velocity dispersion alone. However, we 
have also shown that this degeneracy can be broken by us- 
ing the velocity dispersions in the two schemes simultane- 
ously. In this paper we first tested our method using a mock 
galaxy catalogue. We fitted the measured satellite-weighted 
and host-weighted velocity dispersions simultaneously us- 
ing a parametric model for the halo occupation statistics 
of central and satellite galaxies, and demonstrated that we 
can reliably obtain confidence levels on the true mean and 
scatter of the mass-luminosity relation of central galaxies. 

Next we applied the above method to a volume limited 
sample from the SDSS. The mean of the mass-luminosity 
relation increases as a function of the central host lumi- 
nosity indicating that, as expected, more massive haloes 
host brighter centrals. This result is in excellent quantita - 
tive agreement with a recent study bv lCacciato et al.l ((2003), 



16 More et al. 



who use the abundance and clustering properties of galax- 
ies in SPSS t o con strain the CLF, and with the SAM of 
ICroton et al.1 ()2006l ). The sateUite kinematics obtained in 
our study are consistent with a model in which P{Lc\M) 
has a constant scatter, criogL, independent of the halo mass 
M. We obtain o-iogz, ~ 0.16 ± 0.04 in excellent agreement 
with other independent measurements suggesting that the 
amount of stochasticity in galaxy formation is similar in 
haloes of all m asses. This is also suggested by th e SDSS 
group catal og of Yang. Mo fc van den BoschI (|2008r ) and by 
the SAM of lCroton et al.l (|2006l ). 



ACKNOWLEDGMENTS 

We are grateful to Kris Blindert, Jacqueline Chen, 
Anupreeta More, Peder Norberg, Hans-Walter Rix and 
Joachim Wambsganss for valuable discussion. 



Wang Y., Yang X., Mo H. J., van den Bosch F. C, Wein- 

mann S. M., Chu Y., 2007b, arXiv, 711. larXiv:0711.4431l 
White S. D. M., Rees M. J., 1978, MNRAS, 183, 341 
Yang X., Mo H. J., Jing Y. R, van den Bosch F. C, 2005a, 

MNRAS, 358, 217 
Yang X., Mo H. J., van den Bosch F. C, Weinmann S. M., 

Li C, Jing Y. P., 2005b, MNRAS, 362, 711 
Yang X., Mo H. J., van den Bosch F. C, Pasquali A., Li 

C, Barden M., 2007, ApJ, 671, 153 
Yang X., Mo H. J., van den Bosch F. C, 2008, ApJ, 676, 

248 
York D. G., et al., 2000, AJ, 120, 1579 
Zaritsky D., Smith R., Frenk C, White S. D. M., 1993, 

ApJ, 405, 464 
Zaritsky D., White S. D. M., 1994, ApJ, 435, 599 
Zaritsky D., Smith R., Frenk C, White S. D. M., 1997, 

ApJ, 478, 39 



APPENDIX A: 
GALAXIES 



SAMPLING OF CENTRAL 



REFERENCES 

Adelman-McCarthy J. K., et al., 2006, ApJS, 162, 38 

Becker M. R., et al., 2007, ApJ, 669, 905 

Blanton M. R., et al., 2005, AJ, 129, 2562 

Blumenthal G. R., Faber S. M., Primack J. R., Rees M. J., 

1984, Nat, 311, 517 
Brainerd T. G., Specian M. A., 2003, ApJ, 593, L7 
Cacciato M., van den Bosch F.C., More S., Li R., Mo H.J., 

Yang X., 2008, MNRAS, submitted 
Chen J., 2007a, arXiv, 711. . arXiv:0711. 09891 
Chen J., 2007b, arXiv, 712, larXiv:0712.0003l 
CoUess M., et al., 2001, MNRAS, 328, 1039 
Conroy C, et al., 2005, ApJ, 635, 982 
Conroy C, et al., 2007, ApJ, 654, 153 
Cooray A., 2006, MNRAS, 365, 842 
Croton D. J., et al., 2006, MNRAS, 365, 11 
Davis M., Efstathiou G., Frenk C. S., White S. D. M., 1985, 

ApJ, 292, 371 
Diaferio A., Geller M. J., 1996, ApJ, 467, 19 
Erickson L. K., Gottesman S. T., Hunter J. H., Jr., 1987, 

Nat, 325, 779 

Klypin A., Prada F., 2007, arXiv, 706, arXiv:0706.3554 
Kravtsov A. V., Berlind A. A., Wechsler R. H., Klypin 

A. A., Gottlober S., AUgood B., Primack J. R., 2004, ApJ, 

609, 35 
Maccio A. V., Button A. A., van den Bosch F. C, Moore 

B., Potter D., Stadel J., 2007, MNRAS, 378, 55 
McKay T. A., et al., 2002, ApJ, 571, L85 
More S., van den Bosch F. C, Cacciato M., 2008, MNRAS, 

submitted, arXiv:0807.4529" 
Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ, 490, 

493 
Norberg P., Frenk C. S., Cole S., 2008, MNRAS, 383, 646 
Prada F., et al., 2003, ApJ, 598, 260 

Porciani C, Dekel A., Hoffman Y., 2002, MNRAS, 332, 325 
Spergel D. N., et al., 2007, ApJS, 170, 377 
van den Bosch F. C, Norberg P., Mo H. J., Yang X., 2004, 

MNRAS, 352, 1302 
Wang Y., Yang X., Mo H. J., van den Bosch F. C, 2007a, 

ApJ, 664, 608 



The ultimate goal of satellite kinematics is to probe the halo 
mass-luminosity relationship (MLR) of central galaxies. In 
principle, an unbiased estimate for the MLR requires that 
the central galaxies identified by the selection criterion are 
an unbiased (sub-) sample with respect to their correspond- 
ing dark matter haloes. In this appendix, we investigate, 
using the MGC, how our iterative criterion performs in this 
respect and compare it with the strict criterion used in the 
literature. 

For reasons that will become clear later, we use the 
flux-limited sample MOCKF for this test. The solid lines 
in Fig. I All show the distributions of halo masses, P{M\Lc), 
for all central galaxies in MOCKF divided in 5 luminosity 
bins. The average logarithm of the luminosities of central 
galaxies in each bin is indicated at the top right corner. The 
dotted lines show the distributions, P{M\Lc), for all central 
galaxies that have at least one satellite in the selection aper- 
ture defined by our iterative selection criterion. Finally, the 
histograms show the distributions, P{A4\Lc), of the centrals 
selected by our iterative criterion. Clearly, the centrals se- 
lected by our iterative criterion sample the distribution of 
halo masses from the dotted lines (and not the solid lines). 
However, as discussed in Section [4.31 this bias is taken into 
account while modelling the kinematics (see Eq. [22) ') and 
therefore allows us to make an unbiased estimate. As shown 
in Section |S] we indeed recover an unbiased MLR from the 
kinematics of the selected satellites measured around the 
centrals selected by our iterative criterion. 

For comparison, we now repeat this exercise using the 
strict criterion employed in previous studies. In particular, 
we adopt the criterion used in N08. This criterion identifies 
a galaxy as a central if it is at least /h = 2 times brighter 
than any other galaxy in a fixed (irrespective of the lumi- 
nosity of the galaxy) aperture cylinder (see Table [1]) around 
itself. Satellites are identified as those galaxies that are at 
least /s = 8 times fainter than the centrals and reside in a 
smaller aperture cylinder defined around the centrals. The 
values of /h and /a in the N08 criterion are conservative, as 
the principle goal of their study was to select isolated cen- 
tral galaxies. Applying the N08 criterion to MOCKV selects 
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Figure Al. Comparison of the sampling of central galaxies using the iterative selection criterion (ITER) used in this paper and the 
criterion used in Norberg et al. (2008) (N08). The histograms in the upper (bottom) panel show the distributions of halo masses of 
central galaxies selected according to ITER (N08). The average log(Lc/ /i^'^Lq) for central galaxies in each bin is indicated at the top 
right corner of each panel. The solid lines show the true distributions of halo masses for all the central galaxies and the dotted lines show 
the distribution of halo masses of those central galaxies that have at least one satellite more than /s times fainter than themselves in the 
inner cylinder defined by the selection criterion. 



only 126 satellites around 96 central galaxies. Therefore, to 
do a meaningful comparison, we apply the N08 criterion to 
MOCKF for which it selects 657 satellites around 395 cen- 
trals. For comparison, our iterative criterion yields 39951 
satellites around 21206 centrals. 

Solid lines in the lower panels of Fig. lAll are the same 
as in the upper panels and show the distributions of halo 
masses, P{M\Lc), for all central galaxies in MOCKF di- 
vided in 5 central luminosity bins. The dotted lines show the 
P{M\Lc) for those centrals that have at least one satellite 
around them which is /s(= 8) times fainter than themselves. 
There is a negligibly small difference in the dotted lines in 
the two rows due to different values of /g. Finally, the his- 
tograms show the P(A'I\Lc) distributions of the sample of 
centrals selected by the N08 criterion. Clearly, these do not 
sample the distributions shown by the dotted lines and the 
distributions are clearly biased towards the low mass end, 
especially, in the bright luminosity bins. This owes to the 
fact that Norberg et al. (2008) adopt /h — 2, which pref- 
erentially selects centrals that do not have satellite galaxies 
of comparable brightness. This biases the distributions to- 
wards the low mass end. Note, though, that this is not a cri- 
tique regarding their selection criteria; after all, as Norberg 
et al. (2008) clearly described in their paper, their princi- 
pal goal is to study the kinematics around isolated galaxies. 
However, it does mean that it is not meaningful to compare 
their MLR, which is only applicable to isolated galaxies, to 
that obtained here, which is representative of the entire cen- 
tral galaxy population. 



